Photoproduction of rj and rj ' Mesons off Protons 
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Total and differential cross sections for 77 and 77' photoproduction off the proton have been de- 
termined with the CBELSA/TAPS detector for photon energies between 0.85 and 2.55 GeV. The 
r/ mesons are detected in their two neutral decay modes, 77 — > 77 and r\ — ► 3ir a — > 67, and for the 
first time, cover the full angular range in cos 8 cm of the 77 meson. These new 77 photoproduction 
data are consistent with the earlier CB-ELSA results. The 77' mesons are observed in their neutral 
decay to 7r°7r°77 — > 67 and also extend the coverage in angular range. 

PACS numbers: 11.80.Et,13.30.Eg,13.60.Le,13.75.Gx,14.20.Gk,14.40.Aq,25.20.Lj 



I. INTRODUCTION 

Understanding the structure of the proton and its 
excited states is one of the key questions in hadronic 
physics. Known as the missing-baryon problem, quark 
models based on three constituent quark degrees of free- 
dom predict many more states than have been observed 
experimentally Baryon resonances are broad and widely 
overlap, especially at higher energies, imposing chal- 
lenges on the interpretation of experimental data in terms 
of resonance contributions. Without precise data from 
many decay channels, it will be difficult or even impossi- 
ble to accurately determine the properties of well estab- 
lished resonances, or to confirm or rule out the existence 
of weakly established resonances or new, so-far not ob- 
served states. 

Of particular importance are well-chosen decay chan- 
nels which can help isolate contributions from individual 
excited states and clarify their importance. Photopro- 
duction of i] and i] ' mesons offers the distinct advantage 
of serving as an isospin filter for the spectrum of nucleon 
resonances and thus, simplifies data interpretations and 
theoretical efforts to predict the excited states contribut- 
ing to these reactions. Since the 77 and 77' mesons have 
isospin 1 = 0, the TV 77 and Nrj' final states can only 
originate from intermediate 1=1/2 nucleon states. 

Data on 77 photoproduction off the free proton were 
obtained and studied at many different laboratories over 



a wide kinematic range 0, H !, i, i, i, 0, H ©, 03] ■ A 

review of the main data sets and a corresponding com- 
parison of their coverage in energy and solid angle can 
be found in [ll|. Almost all analyses found that the 
iV(1535)5n nucleon resonance dominates 77 photopro- 
duction at threshold, though there are models which do 
not need the N(1535)Sn resonance to describe thresh- 
old production of 77 mesons [l2j ■ The JV(1535)5n state is 
well-known for its large A 77 coupling, whereas other res- 
onances couple only weakly to Nrj. Small contributions 
from the A(1520)Di3 resonance via interference with the 
Sn resonance have been determined from data on pho- 
ton beam asymmetries 0, and angular distributions. 
Data from target polarization experiments [H revealed 
surprising effects concerning the phase relations of the 
s- and d-wave amplitudes [13| . Despite its four-star as- 
signment by the Particle Data Group (PDG) (HI, the 
role and nature of the A(1535)Sii is still not well under- 
stood. Surprisingly, the AT(1650)5n nucleon resonance 
has the same quantum numbers, but shows no strong TV 77 
coupling. Many different arguments have been discussed 
to explain this observation; the two S*n states can have 
appreciable mixing [l5| . for instance. The AT(1535)5n 
resonance can also be a dynamically-generated state of 
the YjK —prj system fl?jj or more generally, a dynamically 
generated resonance coming from the interaction of the 
octet of pseudoscalar mesons with the ground-state octet 
of baryons [13, EH ■ Recent efforts at Jefferson Labora- 
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tory have concentrated on describing the jp — > N* tran- 
sitions by the interaction of the photon with the 3-quark 
core of the resonance including meson-cloud effects in the 
low Q 2 region The agreement of the model predic- 
tions for the helicity amplitude Ai/ 2 with experimental 
data is good for some lower-lying N* states to fair for 
the -/V(1535)5ii not ruling out alternative explanations. 

The importance of contributions from the N(1650)Sn 
resonance to r/ photoproduction has been discussed fur- 
ther in conjunction with its photoproduction off the neu- 
tron. Recently, the neutron data have attracted inter- 
est due to the observation of a bump-like structure at 
1.67 GeV/c 2 [13, which has not been seen in the 
cross section off the proton. In , it has been shown 
that a strong interference between S l n(1535), 5*11(1650), 
and a non-resonant background can provide a good de- 
scription of these data. 

A partial wave analysis (PWA) of recent CB-ELSA 
data in the framework of the Bonn-Gatchina (BnGa) 
model (23[, which included data on other reactions and 
from other experiments, found the dominance of three 
nucleon resonances in rj photoproduction: ^(lSSS)^! 
iV(1720)Pi3, and a proposed new state, V(2070)Z?i 5 
The large Nrj coupling of the V(1720)Pi3 was surprising. 
Solutions of the ry-MAID model [24[ in this mass range 
assign much of the intensity to the AT(1710)Pn instead. 
Current efforts with regard to the extraction of double- 
polarization observablcs will help shed light on this con- 
troversy. The data presented here cover the full angular 
range and within the framework of the Bonn-Gatchina 
model, are still consistent with the dominance of the 
three nucleon resonances, iV(1535)Sii, V(1720)Pi3, and 
iV(2070)L>i 5 . 

Data on rj' photoproduction is scarce. Analyses pub- 
lished before 2005 observed fewer than 300 r\ ' events [25], 
I2H [27| and an interpretation in terms of resonance con- 
tributions was difficult. Data from Jefferson Laboratory 
significantly improved the world database [28] . They ob- 
served 2 x 10 5 77 ' events, which allowed the extraction of 
differential cross sections. Though more precise than pre- 
vious measurements, the CLAS data are still limited in 
their angular coverage. In the model by Nakayama and 
Haberzettl ([H[ and ref. [12] therein), the Nrj' final state 
couples to AT(1535)5u and N(17W)P n . The authors 
claim the importance of J = 3/2 states (7V(1940)-Pi3, 
N(1780)D 13 , iV(2090)£>i 3 ) in the process, which are use- 
ful to obtain the correct shape of the differential cross 
sections for energies from 1.728 GcV to 1.879 GeV. 

In this paper, we present total and differential cross 
sections for the reactions: 

7P — > prj, where rj — > 27, (1) 
■yp — > pr], where rj — > 3n° — > 67, and (2) 
jp^prj', where rj ' — > 2TT°rj — > 67. (3) 

The data cover an incoming photon energy range up to 
2.55 GeV and show the full angular coverage. 

The paper has the following structure. Section HT1 gives 
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FIG. 1: Experimental setup of CBELS A/TAPS in Bonn. The 
electron beam delivered by the accelerator ELSA enters from 
the left side. 



a brief introduction to the CBELSA/TAPS experimental 
setup. The data reconstruction and selection is discussed 
in section IIIII and the extraction of differential and total 
cross sections is described in section IIVI Experimental 
results are finally presented in section fVl 



II. EXPERIMENTAL SETUP 

The experiment was carried out at the electron acceler- 
ator facility ELSA [2^] at the University of Bonn using a 
combination of the Crystal-Barrel [30[ and TAPS [3l|,|32j| 
detectors. The experimental setup is shown in Fig. [TJ 

Electrons with an energy of 3.175 GeV were ex- 
tracted from ELSA via slow (resonant) extraction. The 
bremsstrahlung-taggcr photon-beam facility at ELSA de- 
livered unpolarizcd tagged-photon beams in the energy 
range from 0.5 to 2.9 GcV by passing the electron 
beam through a thin copper radiator with a thickness 
of (3/1000) • Xr (radiation length). Electrons are de- 
flected in the field of the tagger dipole magnet according 
to their energy loss in the bremsstrahlung process; the 
remaining energy is determined in a tagger detector con- 
sisting of 480 scintillating fibers above 14 scintillation 
counters (tagger bars) in a configuration with adjacent 
paddles partially overlapping. The corresponding energy 



of an emitted photon is E 1 = Eq 



Electrons not 



undergoing bremsstrahlung are deflected at small angles 
and guided into a beam dump located behind the tagger 
detectors. The energy resolution is about 2 MeV for the 
high-energy photons and 25 MeV for the low-energy part 
of the bremsstrahlung spectrum. 

For the energy calibration of the tagger, a polynomial 
was determined in simulations using the measured field 
map of the bending magnet and the known positions of 
the fibers. The calibration was cross-checked by mea- 
surements with the ELSA electron beam at two different 
energies. At 600 and 800 MeV, a low-current beam was 
guided directly into the tagger, while the magnetic field 
was slowly varied. These results provided corrections to 
the initial polynomial [33| . 

The photons hit the liquid hydrogen target in the ccn- 
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FIG. 2: Top: Schematic drawing of the liquid hydrogen 
target, scintillating-fibre detector, Crystal-Barrel and TAPS 
calorimeters. Bottom: Front view of TAPS; the left side shows 
the logical segmentation for the LED-low trigger, the right 
side the logical segmentation for the LED-high trigger (see 
text for more details). 



tcr of the Crystal-Barrel (CB) calorimeter. The target 
cell (5 cm in length, 3 cm in diameter) was surrounded 
by a scintillating-fiber detector (34|, which provided an 
unambiguous impact point for charged particles (due to 
the arrangement of its three layers) leaving the target. 
The CB-calorimeter in its CBELSA/TAPS configuration 
of 2002/2003 consisted of 1290 CsI(Tl) crystals with a 
length of 16 Xfj. The modules have an excellent photon 
detection efficiency; a detailed description can be found 
in pSCj]. For this series of experiments, the (downstream) 
rings 11-13 were removed to combine the detector with 
TAPS in the forward direction. The CB-calorimeter cov- 
ered the complete azimuthal angle and polar angles from 
30° to 168°. All crystals are of trapezoidal shape pointing 
to the center of the target (Fig. [3J top) . 

The TAPS detector consisted of 528 hexagonal BaF2 
crystals with a length of about 12 Xr. It was configured 
as an hexagonal wall serving as the forward end cap of 
the Crystal-Barrel calorimeter (Fig. [3J bottom). TAPS 
provided a high granularity in the forward direction cov- 
ering polar angles between 5° and 30° (full (j> coverage). A 
5 mm thick plastic scintillator in front of each TAPS mod- 
ule allowed the identification of charged particles. The 
combination of the Crystal-Barrel and TAPS calorime- 
ters covered 99 % of the 4ir solid angle and served as an 



excellent setup to detect multi-photon final states. 

The fast response of the TAPS modules provided the 
first-level trigger. The second-level trigger was based on 
a cellular logic (FACE), which determined the number 
of clusters in the barrel. The trigger required either two 
hits above a low-energy threshold in TAPS (LED-low) or 
one hit above a higher-energy threshold in TAPS (LED- 
high) in combination with at least one FACE cluster. The 
shape of the logical segmentation for the TAPS trigger is 
shown in (Fig. [51 bottom). 

The beam-monitor placed at the end of the beam 
line provided valuable information on the beam intensity 
(photons not interacting in the H2 target) used for the 
determination of the photon flux. This total absorption 
Cerenkov counter consisted of an array of 9 lead glass 
crystals. 

III. DATA ANALYSIS 

Data presented here were accumulated from October 
2002 until November 2002 in two run periods with ELSA 
beam energies of 3.175 GeV. These data were used to ex- 
tract differential and total cross sections for a variety of 
final states [35|, HH . The event reconstruction and selec- 
tion of the two r] photoproduction channels ([1]) and ([2]) 
as well as the 77 ' photoproduction channel ([3]) with inci- 
dent photon energies up to 2.55 GeV is presented in this 
section. The total number of ~ 600, 000 77 events was ob- 
served (~ 422, 300 for 77 -> 77 and ~ 126, 300 77 -> 3tt°) 
covering invariant masses from 1510 to 2380 MeV/c 2 . 
For the 77' channel, ~ 5100 events were observed cov- 
ering invariant masses from 1920 to 2380 MeV/c 2 . The 
77' threshold region, M e [1896, 1920] MeV/c 2 , was an- 
alyzed separately using a finer energy binning to better 
study the threshold behavior. 

A. Event Reconstruction 

Events with at most one proton and with two or six 
photons were selected, respectively. The charged clusters 
were identified in TAPS by using the plastic scintillators 
mounted in front of each BaF2 crystal. The efficiencies of 
these (photon)-veto detectors were determined and mod- 
elled in the Monte-Carlo (MC) program. In the Crystal- 
Barrel reconstruction, a cluster is assigned to a charged 
particle if the trajectory from the target center to the bar- 
rel hit forms an angle of less than 20° with a trajectory 
from the target center to a hit in the scintillating fiber 
detector. Proton identification is only used to remove 
it from the list of photon candidates. The proton mo- 
mentum is then reconstructed from event kinematics in 
"missing-proton" kinematic fitting. Proton clusters are 
on average much smaller than photon clusters and pro- 
vide worse angular resolution. The proton momentum 
direction reconstructed from kinematic fitting had to be 
consistent again with a calorimeter hit when a charged 
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Confidence Level (ti -» yy) Confidence Level (11 -» 3ji) Pullofi>(y) Pullof$(y) 



FIG. 3: Distributions of confidence levels resulting from a IC kinematic fit of a true data event sample with two photons in 
the final state to the hypothesis 7p — > P77 (a) and a 4C fit of a 6-photon event sample to the hypothesis 7p — > p37r° (b) 
for £ 7 6 [f600, f650] MeV. The insets show the same for Monte-Carlo data. The strong rise near zero is due to poorly 
reconstructed signal events. Examples of pull-distributions for the fit parameter <fi of the final-state photons are shown in (c) 
for data and in (d) for Monte-Carlo events for the kinematic fit to 7p — > p37r° integrated over the full kinematic range. The 
mean and a values are 0.0036 and 0.00078 as well as a = 1.14 and 1.03, respectively. 



cluster was identified. Our Monte-Carlo studies of reac- 
tion |T]) show that if a proton was observed in the event, 
the overall misidentification probability is less than 3%. 

A kinematic fitter slightly adjusts the measured values 
within the estimated errors by a minimization procedure 
until they fulfill exactly certain constraints expressed in 
the form of equations, which are based on physical condi- 
tions like energy and momentum conservation or invari- 
ant masses. The x 2 probability or confidence level (CL), 
which is derived from the \ 2 value of the fit, can be used 
to make judgements and decisions about the goodness 
of the fit and provides an ideal method to judge possi- 
ble final-state hypotheses for an event. For Gaussian- 
distributed errors of the measured particle properties, 
confidence levels should be flat. Pulls are defined to test 
the correct determination of the covariance matrix and 
are a measure of the displacement of the reconstructed 
values to the fitted values. They are constructed such 
that a valid distribution of pulls will form a normal dis- 
tribution with a width of one and a mean of zero. Pulls 
are very sensitive to the goodness of the fit. If the width 
deviates from one, the resolution derived in the recon- 
struction does not reflect the true errors and it is neces- 
sary to globally scale the measured initial errors in order 



reaction 


decay mode 


constraints 


fit 


7P -> pn 
7P -» pri 

yp->pri' 


r) — > 37T° — > 67 
r] ' — > 27r°?7 — » 67 


(E,p) conservation 
(E.p) conservation 

+ 3 x 7T° mass 
(E,p) conservation 

+ 7T° n 77 masses 


IC 
IC 
4C 
IC 

3C 



TABLE I: Kinematic fits and constraints used in the analysis. 
The proton is treated as missing particle. Its momentum is 
determined from the kinematic fit. For the r\' reconstruction, 
only one pion mass was imposed in the kinematic fit (3C fit). 



to force the pull distributions to have a width of one. 
No scaling factors are needed in this analysis for events 
with two photons in the final state. For events with six 
photons in the final state, scaling factors have been deter- 
mined carefully for data and Monte-Carlo events. Typ- 
ical confidence-level and pull distributions are shown in 
Fig. [3j The CL values were found to be sufficiently flat 
in all photon-energy bins. 

In a first step of kinematic fitting, a consistency check 
was carried out by imposing energy and momentum con- 
servation on all events. The hypothesis 

7P^pn 7 7 (4) 

was tested, where n 7 is the number of photons in the fi- 
nal state, i.e. two for reaction |T]) and six for ^ and 
Energy and momentum conservation provides four equa- 
tions which any event due to reaction ([4]) must satisfy. 
Thus, the proton three-momentum can be left "missing" 
and reconstructed from other observables, still retaining 
one constraint (IC) provided by the photon energies and 

xio! 
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FIG. 4: Time difference between photons in TAPS (mean 
value) and electrons in the tagger integrated over all events. 
A prompt coincidence is defined by —3 < At < 3 ns. 
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FIG. 5: Invariant mass spectra (integrated over all photon energies) for the reactions 77} — » P77 (left) and jp — * p-K°n°iv 
(center); confidence-level cuts were applied at 10 -2 and at 10 -3 , respectively. In the 2-photon decay mode, the 7r° and r\ 
mesons are observed with very little background. The invariant 7r°7r°77 mass spectrum (right) shows a clear r\' signal and an 
enhancement at 1250 MeV/c 2 . 



directions. Tabic Q] summarizes the hypotheses used to 
select events for jp — > prj and jp — * prj 1 . 

A prompt coincidence between a photon in TAPS and 
an electron in the tagger was required to reduce time 
accidental background. Random time coincidences un- 
derneath the prompt peak (Fig. 2]) were subtracted by 
performing the exact same selection procedure for events 
outside the prompt time coincidence window. Fig. [5] 
(left) shows the invariant 77 mass (time-accidental back- 
ground subtracted) for kincmatically-fittcd two-photon 
events. A confidence level cut at 10 ~ 2 was applied. 
Clear peaks for the ir° and 77 mesons are visible. The 
background underneath the peaks depends on kinemat- 
ics and is on the average smaller than 4 %, but can be up 
to 15% at high energies and in the most forward angle 
bins. In addition to energy and momentum conserva- 
tion, three it mass constraints were imposed on events 
due to 7p — > prj — > p37r° — > p6j (4C kinematic fit) reduc- 
ing significantly combinatorial background. A confidence 
level cut at 10 -3 was applied. The invariant 3ir° mass 
is shown in Fig. [5] (middle). The 77 peak is visible above 
a small combinatorial background. The uncertainty in 
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FIG. 6: Invariant 77 mass spectra (integrated over all incom- 
ing photon energies) for the reactions 7P — » p7r°7j77 (left) and 
7P — > p7r°7r°77 (right). The ty° signal shows a double-peak 
structure since the "better" tv° was found in the kinematic fit. 



signal lost due to CL-cuts in data compared to Monte 
Carlo is estimated to be less than 3 %. 

In a final step of the 77 analysis, mass cuts were applied 
in the two-photon spectrum and in the 37r° spectrum. 
The width of the 77 peak varied as a function of incident 
photon energy between ~ 12 MeV/c 2 at the lowest ener- 
gies and ~ 17 MeV/c 2 at the highest energies. The re- 
maining background underneath the 77 77 and 77 37r o peaks 
was subtracted using side bins for every (E 1 , cos 8) bin. 

For the selection of 77' events, the hypothesis ^p — > 
p7r°?777 was tested in addition to reaction with a 
confidence-level cut at 10~ 2 . The remaining invariant 
77 mass is shown in Fig. [6] (left). The second ir° was 
reconstructed with 110 < m jr o_ >Tr < 160 MeV/c 2 . Fig. [5] 
(right) shows a clear peak for the 77 ' in the invariant 
27r°r7 mass spectrum. Moreover, an interesting enhance- 
ment is visible at 1250 MeV/c 2 giving rise to a possible 
contribution of the controversial meson Ty(1295) and/or 
the /i(1285). Further studies of this signal are statis- 
tically challenging and are not discussed further here. 
Events due to reaction ([3]) were finally selected with 
910 < m 27r o^ < 1010 MeV/c 2 . The remaining back- 
ground was determined in fits to the 77 ' peaks. An alter- 
native way of reconstructing 77 ' events via 775 — > p7r°7r°77 
was used for systematic checks. The invariant 77 mass 
is shown in Fig. O (right). 



B. Monte-Carlo simulations 

The performance of the detector was simulated in 
GEANT3-based Monte-Carlo studies. The program 
package used for CBELS A/TAPS is built upon a program 
developed for the CB-ELSA experiment. The Monte- 
Carlo program reproduces accurately the response of the 
TAPS and Crystal-Barrel crystals when hit by a photon. 
For charged particles, the detector response is known to 
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a lower precision but still reasonably well understood. 

The acceptance for the reactions (H}-© was deter- 
mined by simulating events, which were evenly dis- 
tributed over the available phase space. The Monte-Carlo 
events were analyzed using the same reconstruction crite- 
ria, which were also applied to the (real) measured data. 
The same hypotheses were tested in the kinematic fits 
and events selected with the same confidence level cuts. 
The acceptance is defined as the ratio of the number of 
generated to reconstructed Monte-Carlo events 

A lv _ vX = -r T (X = T), j] ). (5) 

Jv gen,MC 

The trigger required either two hits above a low-energy 
threshold in TAPS (LED low - Leading-Edge Discrimina- 
tor threshold) or one hit above a higher-energy threshold 
in TAPS (LED high) in combination with at least one 
cluster in the Crystal-Barrel calorimeter. The second- 
level trigger used a fast cluster encoder (FACE) based 
on a cellular logic to define the number of contiguous 
Crystal-Barrel clusters. The decision time depended on 
the complexity of the hit distribution in the Crystal- 
Barrel and was typically 4 [is. In case of an event rejec- 
tion, a fast reset was generated, which cleared the read- 
out electronics in 5 [is. Otherwise the readout of the full 
event was initiated with typical readout times of 5-10 ms. 
To properly simulate the detector response, FACE and 
TAPS-LED thresholds had to be determined from the 
data for all crystals. Given the different response charac- 
teristics of protons and photons in BaF2 crystals, protons 
experience slightly different LED-thrcsholds than pho- 
tons. For this reason, we have corrected the measured 
proton energy according to 0.8 • E p + 30 McV, which is 
derived from available proton times in TAPS and from 
Monte-Carlo studies. At the reaction threshold, when 
the proton is required in the (TAPS) trigger, corrections 
are small. Above about 1 GeV in the incoming photon 
energy, the proton trigger is not relevant. Our under- 
standing of the threshold function is fair to good and 
reasonably well reproduced in the trigger simulation. 



The target area density, i.e. the number of atoms in 
the target material per cross-sectional area (orthogonal 
to the photon beam), is given by 

^ = 2 # e1 7tTT- 2 - 231 - 10 ^ 1 ' W 

M mo i(t± 2 ) 

where p(H 2 ) = 0.0708 g/cm 3 is the density and M mo [ = 
2.01588 g/mol the molar mass of liquid H 2 . N A = 6.022 • 
10 23 mol -1 is the Avogadro number and L = 5.275 cm 
the length of the target cell. The factor of two accounts 
for the molecular composition of hydrogen (H2 ) . 

The cross sections were extracted independently for 
both r\ decay modes, r\ — > 77 as well as r\ — > 37r°, and 
then averaged (weighted with their errors) based on the 
observed good agreement (Fig.[7])- The number of events 
in an (£7, cos 9 cm ) bin comprises events with two or three 
final-state particles (at least 27's) or six or seven "par- 
ticles" (at least 67's), respectively. The proton can be 
"missing" , but events with and without detected proton 
are treated in the same way in the event reconstruction. 
At threshold, the event kinematics requires that the pro- 
ton is used in the (TAPS) trigger. Thus, the threshold 
function for the detection of low-energy protons needs to 
be reasonably well understood. 

The interval of the solid angle is given by Af2 = 
27rA cos(0 cm ) with A cos(# cm ) describing the width of the 
angular bins. It was chosen to be 0.1 for the r\ data 
presented here. Energy bins were defined by consider- 
ing statistics and ensuring a good comparability with 
other experiments. A total of 34 bins is presented in 
energy steps of 50 McV for E 1 E [850, 2550] MeV. For 
the rj' data, Acos(£? cm ) was chosen to be 0.2. A total 
of 20 bins is presented in energy steps of 50 MeV for 
E 7 e [1500, 2550] MeV. 

The number of observed rj and 77' mesons needs to be 
corrected for unseen decay modes. Partial-decay branch- 
ing fractions used to correct the measured cross-sections 
arc taken from Q: BR(?j -> 2 7 ) = 0.3931 ± 0.002, 
BR(r? -> 3tt° -> 67) = 0.3256 ± 0.0023, and BR(?/' -> 
27r°7/ -> 67) = 0.207 ±0.012. 



IV. DETERMINATION OF CROSS SECTIONS 

The differential cross sections for this analysis are de- 
termined according to 

da _ JV^„ l7 1 1 1 

dQ A x ^ n ^ N lPt Ail ££^r' 1 > 

Ttotal 

where 

pt : target area density 
Nx^„ 7l : number of reconstructed data events 
in an (i5 7 , cos # cm ) bin 
iV 7 : number of photons in an (-E 7 ) bin 
AA-^n T 7 : acceptance in an (iJ 7 ,cos# cm ) bin 

A£7 : solid-angle interval AO = 2itAcos (6 cm ) 

p^" t7 : decay branching fraction. 



A. Normalization 

The tagging hodoscope consisted of 480 scintillating 
fibers above 14 partially overlapping scintillation coun- 
ters (tagger bars). The photon flux was measured di- 
rectly in the experiment and determined according to 

N 7 = tf£Er-a-P 7> (8) 

where N are the free hardware counts for the indi- 
vidual fibers corrected for the life-time of the detector. 
The parameter a accounts for the (fiber)-cluster recon- 
struction in the tagger, which has to be performed in the 
same way as for real hadronic events. The photon def- 
inition probability or P 7 denotes the probability that a 
real photon is emitted along the beam axis in the tagger 
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for each (_E 7 , cos # cm ) bin. 



The values are weighted with their squared inverse errors. 

and traverses the liquid hydrogen target. The scalers are 
recorded in scaler events, which were accumulated with a 
minimum-bias trigger at a rate of 1 Hz during the regular 
data taking. This trigger required only a hit in the tag- 
ger and was thus independent of hadronic cross sections. 
The parameter P 1 is determined from Tagger- Or- Runs — 
separate data runs utilizing a minimum-bias trigger. 

The total error of the photon flux is assumed to be 
dominated by the P 1 error and depends strongly on the 
efficiency of the beam monitor at the end of the beam line 
(Fig. H]). P 7 was determined to 0.639±0.002 stat ±0.05 sys . 
by varying the background subtraction of non-coincident 
tagger-beam monitor hits. This value is consistent with 
determinations from multiple Tagger- Or- Runs at differ- 
ent incoming photon rates. An overall error of 10 % has 
been assigned to the photon flux determination. 



B. Systematic Uncertainties 

The statistical errors are determined from the number 
of events in each (E 7 ,cos9 cm ) bin. Statistical errors are 
shown for all data points; systematic uncertainties are 
given as error bands at the bottom of each plot. 

Sources of the systematic errors are uncertainties in the 
exact position of the liquid hydrogen target and a possi- 
ble offset of the photon beam. The position of the target 
cell was determined from kinematic fitting by comparing 
the off-zero displacement of different pull distributions to 
Monte-Carlo simulations. It was found to be shifted up- 
stream by 0.65 cm (37l | . The corresponding systematic er- 
rors were determined by varying the target position in the 
Monte Carlo (±1.5 mm) and evaluating changes in the re- 
extracted differential cross sections. The errors show an 
angular dependence, but are 2-3 % on the average and 
< 5 % at most around cos 9 cm = 0. The photon beam 
was assumed to be shifted by less than 2 mm off axis at 
the target position. The errors of the decay branching 
fractions are negligible. The uncertainty of the proton 
trigger has been determined from the small disagreement 
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FIG. 8: Differential cross sections for jp — » pn using the 
n — > 27 (■) and 77 — > 3ir a — > 67 (a) decay channels. Only 
statistical errors are assigned to the data points. In some 
bins, the acceptance for n — > 37r is small (< 5%) and thus, 
corresponding data points are not shown. The figure shows 
the excellent agreement of both reactions. 



of the differential rj cross sections using the rj —> 2j and 
r/ — > 37T° — > 67 decay channels for E 1 < 1 GeV and 
cos 6 cm < 0.0 (Fig.©. 

The reconstruction of neutral mesons and the identifi- 
cation of final states requires a sequence of cuts including 
the use of kinematic fitting. As discussed in the following 
section, the reconstruction of rj mesons from final states 
with two and six photons leads to compatible results. 
This fact emphasizes a good understanding of the de- 
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FIG. 9: Differential cross sections for 7p — > pr\ (■) using the combined data set of r\ — > 27 and r) — » Sir — * 67 events. 
For comparison, CB-ELSA data || are represented by (a) and CLAS data by (*)■ The solid line shows our previous PWA 
solution 0| and the dashed line represents the SAID model |3^|. The data points include statistical errors only; the total 
systematic error is given as error bands at the bottom of each plot. 



tector response to multi-photon final states. An overall 
±5.7% error is a ssig ned to the reconstruction efficiency 
as determined in [381 ] . An additional 3 % systematic error 
accounts for the slightly different effects of confidence- 
level cuts on data and Monte-Carlo events. All these 
errors are added quadratically to give the total system- 
atic error. Moreover, the 77' systematic error receives an 
additional contribution from an alternative way of recon- 
structing events via jp — > p7r°7r°77 (Sec. IHI Al) . 



V. EXPERIMENTAL RESULTS 

A. Differential Cross Sections dcr/dtt for 7p — > pr? 
at an Electron-Beam Energy of E e - = 3.18 GeV 

Fig. [5] shows the 7p —> prj differential cross sections 
for the two different 77 decay modes ([T]) and ([2]). We 
have excluded those data points in the analysis showing 
a Monte-Carlo acceptance of < 5%. The data sets show 
excellent agreement. 
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dashed line represents the SAID model [3!|. 



We have checked the consistency of the two cross- 
section measurements by forming the ratio of partial 
widths of the two rj decay modes. Since the cross sections 
in Fig. [8] are corrected for the decay branching ratios, we 
determine 





j^(-E 7 ,cos 6 cnl ) 


Tj— »37T° 


r 


- PDG 






^(.E 7 ,cos0 clIl ) 






PDG 



(9) 



for each (-E 7 , cos# cm ) bin. The values are weighted with 
their squared inverse errors and histogrammed (Fig. \T§. 
We derive a peak position of 

^ 3 "° = (0.8133 ±0.0006 sta t. ± 0.0138 sys .) , (10) 
1 1)^27 

where the systematic error is derived from considering all 
systematic uncertainties discussed in Section IIVBI We 
do not claim a new measurement of 77 branching fractions 
here, but have rather used this number to check our re- 
construction efficiency. The Particle Data Group gives 
two values for this branching ratio [lj]: (0.829 ±0.007) is 
the mean value of all direct measurements while a com- 
bined fit to all partial decay widths yields a value of 
(0.828 ± 0.006). Both PDG values are consistent with 
our result. 

B. Differential Cross Sections drr/dfl for 7p — > pr\ 
— Combined Data Set — 

Since the differential cross sections for the two 77 decay 
modes are consistent, we have calculated error- weighted 



mean values. These are presented in Fig. [9] as functions 
of energy and the 77 production angle. Since the cross 
sections change rather smoothly, only few resonances are 
likely to contribute to the process. The N(1535)Sn state 
is known to dominate the threshold region resulting in a 
flat distribution in cos# cm ; interference with other ampli- 
tudes leads to deviations from that flat distribution. At 
higher energies, above E 1 = 1.5 GeV, the development 
of a forward peak indicates important contributions from 
t-channel p and uj exchange. 

Fig. [5] also shows a comparison of our new 77 results 
to published results from CB-ELSA @- The agreement 
between the two data sets is very good at lower ener- 
gies. However, the differential cross sections reported by 
CB-ELSA show somewhat larger discrepancies at higher 
energies and forward angles. Above E 1 = 2.5 GeV and 
in the forward most angle bins, CB-ELSA results are ap- 
proximately 30% larger than our new CBELS A/TAPS 
results at these energies, but are still consistent within 
the errors. We believe that the discrepancy is due to 
underestimated background in the CB-ELSA data in 
the low-statistics forward most bins where fewer than 
10 events were observed. 



C. Differential Cross Sections dcr / dfl for 7p — > prj ' 
at an Electron-Beam Energy of E e - = 3.18 GeV 

Fig. [UJ1 shows the differential cross sections for the re- 
action jp — > prj 1 — > 27r°?7 (reaction [3]) using cos 6£ m< 
bins of width 0.2. The data cover the full angular range. 
Very similar to 77 photoproduction, a rather flat angu- 
lar distribution is observed at low energies suggesting 
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FIG. 11: Differential cross sections for the reaction 7p — -> prj' — > 2n a rj — > p67 close to the reaction threshold (■) determined 
for individual tagger channels using 5 cos 9£. m . bins of width 0.4. Energies in the plots are given in MeV. Though limited in 
statistics, all angular distributions appear to be flat indicating s-wave behavior of the reaction at the threshold. 



s-channel resonance production near threshold. Both 
data sets also show a continuing increase in slope at for- 
ward angles, which becomes more prominent at higher 
energies. This forward peak is most likely due to t- 
channel exchange mechanisms. Moreover, our new data 
indicate a decrease in the forward most bin, which has 
not been observed before. Above 2 GeV in photon en- 
ergy, growth at backward angles could be indicative of 
u-channcl contributions. The overall agreement between 
the CBELS A/TAPS and CLAS data is good at threshold 
to fair above = 1800 MeV. 

The s-wave behavior of the reaction close to the 
reaction threshold is apparent from the experimental 
data. Fig. [IT] shows the differential cross sections for 
'yp — ► pij 1 at and close to the reaction threshold using 
5 cos0 c ' , m bins of width 0.4. The cross sections have 
been determined for individual fibers of the tagging sys- 
tem and cover the full angular range. The data points 
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FIG. 12: Shown is the linear energy dependence of the squared 
total cross section, af ot , for the reaction ~/p — > pr]' close to 



suffer from low statistics, but are overall consistent with 
flat angular distributions. The expected energy depen- 
dence of the reaction at the threshold is given by [40j : 



<7{E 1 )^{E 1 ~E thrcs ) l+1 ' 2 1 = 0,1,2, 



(11) 



the reaction threshold of E*- 



1447 MeV. 



where I denotes s, p, d... waves, etc. For s-wave dom- 
inance, a linear energy dependence of the squared total 
cross section is thus expected close to the reaction thresh- 
old. The differential cross sections shown in Fig. I 111 have 
been used to determine the total ij ' cross section and to 
study the energy dependence. Since the data cover the 
full angular range, no extrapolation is needed. Fig. [T^] 
shows the incoming photon energy plotted versus the 
squared total cross section. A linear dependence is clearly 
observed and a fit determines the energy threshold to 
(1442.6 ± 3.8) MeV, which is compatible with the value 
of (1446.38 ± 0.48) MeV derived from the rj' mass listed 
in the PDG [l4j]. The mass of the rj meson was deter- 
mined by the TAPS/A2 Collaboration in a very similar 
procedure [i"l.l4ll]. 



D. The total cross section 

Fig. [T3] shows the total cross section for r] photopro- 
duction. Due to the complete solid angle coverage, no 
extrapolation is required and the data points are truly ex- 
perimental. In the low-energy range, the 5*11 partial wave 
dominates the cross section. The solid line represents our 
previous PWA solution and is not a fit to these data; the 
two states V(1720)Pi 3 and AT(2070)£>i 5 saturated the 
total cross section Q. A new PWA solution including 
the new CBELS A/TAPS data presented here and other 
data sets is in preparation. It is clear that single- and 
double-polarization variables are required to firmly es- 
tablish resonance contributions. Coupled channel fits to 
many reactions can also help; in particular, when three- 
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body final states arc included, the phase of two-particle 
partial-wave amplitudes is tested in the crossed channel. 

A small anomaly is observed at 1.73 GeV/c 2 in the 
total rj cross section. As discussed earlier, recent data 
off the neutron show a pronounced bump-like struc- 
ture at 1.68 GeV [21j . which has been suggested to 
signal the existence of a narrow baryon state with 
(M « 1.68, T< 30)MeV/c 2 . In particular, the possibil- 
ity that this state could be the iV(1650)Pn, non-strange 
member of an anti-decuplet of pentaquarks is certainly 
interesting [H, H3, H3]. Fig. [14] shows the total number 
of rj — > 77 events for the tagger channels 419-431. The 
data have been fitted using a polynomial to indicate the 
smooth behavior of the distribution. The data point at 
1.73 GeV/c 2 in the total cross section is based on the pho- 
ton energy interval E 1 G [1100, 1150] MeV defined by 
the tagger channels E~ t € [421,426]. No statistically sig- 
nificant enhancement is observed in Fig. ll4l over the small 
energy range under investigation to explain the anomaly 
and a narrow state compatible with the observation in 
the total cross section can be ruled out. We believe that 
this structure is an instrumental effect originating from 
tagger cluster-size corrections of low-rate Tagger-Or-data 
used in the photon-flux determination, which occurs only 
for those channels. 

In Fig. [TS] the total cross section for rj' photopro- 
duction is displayed. Again, due to the complete solid 
angle coverage, no extrapolation is required. At 2 GeV 
in the invariant mass and above, the cross section for 
7?' production is about 50% of that for rj production 
assuming that the processes are dominated by p, u- 
exchange resonances. For a pscudoscalar mixing angle 
of Ops = —19.3°, the non-ss components of r\' and r\ 
differ by a factor 1/2. The similarity of the two numbers 
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FIG. 13: Total 7p — > prj cross section. The data points (•) 
are calculated by summation of the differential cross section, 
the grey line represents the result of our previous partial wave 
analysis. 
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FIG. 14: Total rj — » 77 yields per tagger channel number cov- 
ering the energy range of the anomaly observed in the total 
rj cross section. The data point at 1.73 GeV/c 2 (see Fig. I13[) 
is based on channels 421-426. No statistically significant en- 
hancement can be seen in the excitation function. 



suggests that the dynamics of rj and rj ' photoproduction 
is similar. 



VI. RECENT RESULTS OF THE 
BONN-GATCHINA MODEL 

The analysis of previous CB-ELSA data on photopro- 
duction of 77-mesons [ll| revealed two surprises: 

(1) The analysis suggested a new resonance with spin 
and parity J p = 5/2", 7V(2070)L»i 5 [13]. The total cross 
section was nearly saturated with three resonances, the 
well-known N(1535)Su, the AT(1720)Pi 3 , and the new 
resonance. 

(2) The strong coupling of N (1720)Pi 3 -> Nrj was also 
unexpected; in MAID 0, [H, the iV(1710)P u -> Nrj 
was very significant while A^(1720)Pi3 hardly contributed 
to rj photoproduction. 
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FIG. 15: Total 7p — > prj' cross section. The data points (•) 
are calculated by summation of the differential cross section. 



12 



The pattern of states contributing most to 77 photo- 
production, iV(1535)S u , 7V(1720)Pi3, 7V(2070)D 15 , was 
interpreted as a sequence of quark model states with total 
intrinsic angular momenta L = 1,2,3 and 5 = 1/2 cou- 
pling to J p = 1/2", 3/2", 5/2". The regularity of this 
pattern was used to argue that the seed of AT(1535)5n 
should be of a three-quark nature; due to the presence of 
S- wave thresholds, the state may attract large Nrj and 
E-7T molecular components. 

The relative strength of the two nucleon excitations 
at an incoming photon energy of 1700 MeV in ^yp — > prj 
remains disputed. The GieBen group [46| found - like 
MAID - iV(1710)Pn — > Nrj to provide a significant con- 
tribution. Other coupled channel analyses confirmed 
the dominance of jV(1720)Pi 3 (relative to A(1710)P n ). 
In [47J and [48[, a large variety of r] production data 
was fitted using an effective Lagrangian approach; in 
both analyses, the iV(1720)Pi3 contribution was con- 
siderably larger than that of AT(1710)Pn. A chiral 
quark model approach complemented with a one-gluon 
exchange model [4£| arrived at the same conclusion. 

Restricted to photoproduction data, the best solution 
for the new data presented here (in terms of its \ 2 value) 
still supports the dominance of the three nucleon reso- 
nances, A(1535)5n, A(1720)Pi 3 , and A(2070)£>i 5 in T) 
photoproduction. However, this solution is incompatible 
with data on ir~p — > nr]. This is presently investigated 
further and will be subject of a forthcoming publication 
of the Bonn-Gatchina partial wave analysis group. 



VII. SUMMARY 

In summary, we have presented data on the photo- 
produced r] and r] ' cross sections from the reactions jp — > 
prj with r\ — > 3ir° — > 67 as well as 77 — > 27 and from 
the reaction ^p — > prj' with T} 1 — > 7r°7r°r/ — > 67. The 
continuous beam from the ELSA accelerator and the fiber 
detector of the tagging system provided tagged-photons 
in the energy range from 850 to 2550 MeV. The results 
are in very good agreement with previous measurements, 
but extend over the full angular range in cos # cm of the 
77 and 77 ' meson. The inclusion of the new 77 data into 
a multi-channel partial wave analysis is in preparation. 
The threshold behavior of the 77 ' data indicate s-wave 
dominance. 
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